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ABSTRACT 

New statistical properties of dark matter halos in Lagrangian space are pre- 
sented. Tracing back the dark matter particles constituting bound halos resolved 
in a series of N-body simulations, we measure quantitatively the correlations of 
the proto-halo's inertia tensors with the local tidal tensors and investigate how 
the correlation strength depends on the proto-halo's sphericity, local density and 
filtering scale. It is shown that the majority of the proto-halos exhibit strong 
correlations between the two tensors provided that the tidal field is smoothed 
on the proto-halo's mass scale. The correlation strength is found to increase as 
the proto-halo's sphericity increases, as the proto-halo's mass increases, and as 
the local density becomes close to the critical value, 5ec- It is also found that 
those peculiar proto-halos which exhibit exceptionally weak correlations between 
the two tensors tend to acquire higher specific angular momentum in Eulerian 
space, which is consistent with the linear tidal torque theory. In the light of our 
results, it is intriguing to speculate a hypothesis that the low surface brightness 
galaxies observed at present epoch correspond to the peculiar proto-halos with 
extreme low-sphericity whose inertia tensors are weakly correlated with the local 
tidal tensors. 

Subject headings: cosmology: theory — large-scale structure of universe 

1. INTRODUCTION 

Galaxies are biased tracers of the underlying dark matter distribution. It has been well 
known that the two-point correlation function of the observed galaxies follows a power-law. 
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differ ent from that of the dark matter determined in N-body simulations (e.g.. lMaddox et al. 
1990l ). To test theoretical predictions based on the dark matter against real observations 
of galaxies, it is required to determine a hidden connection between the galaxies and dark 
matter. The density peak formalism was the first attempt to provide such a connection, 
according to which the galaxies form in the high-peaks (i.e., local rnaxirna) of the initial 
smoo thed density field (jPavis et al.l Il985l : iBardeen et al.l Il986l : iKaiserl Il986l : iBond fc Myers 
19961 ). Fitting quite well into the cold dark matter paradigm and explaining successfully some 
observed properties of the galaxies, this density peak formalism became the most prevalent 
model for biased galaxy formation. Given that gravity is mainly responsible for the formation 
and evolution of the galaxies, it was indeed reasonable, appropriate, and natural to regard 
the initial density peaks as good indicators of the sites of galaxy formation. 



An uncomfortable truth, however, had to be faced when iKatz et al.l (119931 ) reported the 
unexpected results derived from N-body simulations that the dark matter particles of galactic 
halos are not well overlapped with those from the high-peaks of the initial density field. They 
traced the trajectories of dark matter particles that comprise the initial density peaks and 
found that the particles from the initial peaks do not end up in real halos. Furthermore, it 
was also shown that the correlation function of high density peaks is different from that of 
galactic halos. These disturbing results forced them to conclude that the high-peaks of the 
initial density field cannot be good indicators of the sites of galaxy formation. In the same 
spirit, iPorciani et al.l (l2002bl ) have shown that only ~ 40% of the proto-halos in a N-body 
simulation contain a density peak within their Lagrangian volume and that gravitational 
shear plays an important role in shaping the proto-halos. 

To find a good indicator of the initial sites of galaxy formation, it is first necessary 
to understand the statistical properties that the majority of the proto-ga lactic sites p ossess 
in Lagrangian space. Such a property was first found unwittingly by iLee fc Peril (120001 ) 
while studying the origin of the galaxy angula r momentum. According to the linear tidal 
torque theory (IDoroshkevichl 1 19 70l : IWhitelll984l ). the angular momentum of a proto-galaxy is 
generated at first order only when the local tidal shear tensor is not perfectly correlated with 
the inertia momentum te nsor of the proto-galaxy. In previous works dealing with the linear 



tidal torque theory (e.g .. ICatelan fc Theunslll996l ). it was assumed that the two tensors are 



generally uncorrected. iLee fc Peril (120001 ) tested the validity of this assumption against the 
numerical data from N-body simulations. They calculated the correlations between the two 
tensors at proto-halo sites and found for the first time that the two tensors are in fact quite 
strongly correlated with each other in contrast to the general assumption. Their result, 
however, was based on rather low-resolution N-body simulations and thus failed to draw 
serious attentions. 
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Later, iPorciani et al.l (l2002bl ) re-tested this assumptio n against hi K h-reso lution simula- 
tions, confirming and extending the prehminary results of iLee fc Peru ( 120001 ) . Noting that 
the existence of strong correlations between the two tensors is such a comrn on phenomenon 
exhibited by most of the proto-halos they considered, iPorciani et al.l fl2002bl ) have suggested 
that the boundaries of the proto-halos are determined by the push and pull of the external 
mass distribution. This is in some sense the opposite of the common wisdom where the key 
factor is assumed to be the self-gravity attraction. 



In this paper we want to further improve upon the previous work by iLee fc Peru (120001 ) 
and IPorciani et al.l (l2002bl ) by using a more detailed statistical treatment and simulations of 
better quality. Our goal here is to quantify the correlations between the two tensors as new 
Lagrangian statistics of dark halos and investigate how the correlation strengths depend on 
proto-halo's shapes and local environments. The outline of this paper is as follows. In §2, 
a brief description of N-body data is provided and the statistical analysis of it is presented. 
In §3 is presented new Lagrangian statistics of dark halos related to the correlations of the 
proto-halo's inertia tensors with the local tidal tensors. In §4 an implication of our result on 
the low surface brightness galaxies is explained. In §5 the results are discussed and a final 
conclusion is drawn. Throughout this paper, we assume a spatially fiat ACDM cosmology. 



2. DATA AND ANALYSIS 



We use the sam ples of dark r natter halos obtained from three different N-body simu- 
lations conducted by iHahn et al.l (120071 ) in periodic boxes of linear size Li = 45/i~^Mpc, 
L2 = 90/i~^Mpc and L3 = 180/i"^Mpc. All three simulations assume a spatially flat ACDM 
cosmology with Q^n = 0.25, fl^ = 0.75, fib = 0.045, as = 0.9, Hq = 0.73 and Ug = 1, each 
following the evolution of 512^ particles to the present epoch from a given initial redshift Zi 
(Li, L2 and L3 have Zi = 79, Zi = 65 and Zi = 52, respectively). Bound halos were resolved 
in each simulation by applying the standard friends-of-friends algori thm with linking-lengt h 
parameter ofb = 0.2/ here / denotes the mean inter-particle distance (lEfstathiou et al.lll985l ). 
Among them, only those halos comprising more than 300 dark matter particles were selected 
to avoid possible numerical artifacts. A total of 50839 halos are selected (13390, 16339 and 
21110 halos from the Li, L2 and L3 simulations, respectively), which span a wide mass range 
of [10.2, 15.3] in units of h~^MQ. A full descr i ption of the simulations and the process of the 
halo-identification is provided in iHahn et al.l (120071 ) . 



Tracing back to the initial conditions the trajectories of dark matter particles that consti- 
tute each selected halo, we locate the proto-halo sites in the Lagrangian space corresponding 
to the initial epoch Zi. Then, we determine the center of mass of each proto-halo site using 
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the positions of its constituent particles. At each halo's center of mass, we measure the mean 
density contrast 6 in the initial density field smoothed with a top-hat filter of scale radius 
Rs- Here we consider four different filtering scales: Rs = 0.5, 1, 2 and 5 h~^Mpc. The initial 
peculiar potential field was derived from the density field by solving the Poisson equation 



( iHahn et al.l 120071 ). and the tidal shear tensor T = (Tij) was measured as the second deriva- 
tive of at the center of mass of each proto-halo. Diagonalizing T, we determine the three 
eigenvectors {ti,t2,t3} corresponding to the three eigenvalues {Ai,A2,A3} in a decreasing 
order. 

The inertia momentum tensor I = (lij) of each proto-halo site is also deter mined using 



the po sitions of the constituent particles in accordance with the formula given in lHahn et al. 
fl2007[ ): 

a 

where m^, is the mass of a-th particle, is the position vector of the a-th particle from the 
center of mass of the proto-halo site, and 6jk is the Kronecker symbol. When the correlations 
with the tidal field smoothed on the scale Rg are calculated, we consider only those halos 
whose mass M is in the range of 0.9Ms < M < I.IM^ where Ms is the top-hat mass enclosed 
by the radius Rs, given the fact that a given proto-halo with mass M would exp erience the 



stron gest effect from the tidal field smoot h ed on the same mass scale M = Mg ( iLee fc Pen 



2OOOI . see also the appendix in iHahn et al.l (120091 )). Table [T] lists the number (Nh) and mean 



mass (M) of halos considered for the correlations with the tidal field smoothed on the four 
different filtering radii Rs- 

A similarity transformation is performed to reexpress I in the principal frame of T. To 
quantify the degree of the correlations between the two tensors, we define a parameter /3, as 

where {qu, Q22, Q33} and {qu, Q23, Q31} represent the three diagonal and off-diagonal elements 
of I in the principal frame of T, respectively. If a proto-halo region has a perfectly spherical 
shape, then the eigenvectors of its inertia tensor are degenerate (i.e., any axis frame can 
be its eigenvector system) and all of the off-diagonal elements of I are always zero (^12 = 
Q23 = ^'31 = 0) .Thus for the case of a perfectly spherical proto-halo region, we always 
have /3 = 1. When a proto-halo region is not perfectly spherical but its inertia tensor I 
is perfectly correlated with the tidal shear tensor T measured at its center of mass, then I 
should be completely diagonal in the principal axis frame of T. Thus for the case of a perfect 
correlation between I and T, we will also have (3=1. On the other hand, if the two tensors 
are uncorrelated, then I in the principal axis frame of T is not diagonal and the off-diagonal 
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elements should be as large as the diagonal ones in magnitude unless the eigenvectors of I 
are degenerate. As the strength of the correlations between I and T decreases, the degree of 
the deviation of jS from the value of unity will increase. 

In the following section, we determine the probability distribution of (3 and investigate 
how the value of (3 depends on the proto-halo's shape, local density and filtering radius. 



CORRELATIONS BETWEEN INERTIA AND TIDAL TENSORS 



Using the numerical data described in §2, we first determine the probability density 
distribution, p{P). Figure [1] shows for the four different cases of Rs = 0.5, 1, 2 and 
5 /i~^Mpc (top-left, top-right, bottom-left, and bottom-right panel, respectively). It is worth 
mentioning again here that when we calculate (3 based on the tidal field smoothed on the 
scale Rg we consider only those proto-halos whose masses belong to the range (O.QMg, I.IM.,) 
with Ms = {A'K/?,)pRl, where p is the mean mass density of the universe. In each panel the 
errorbars represent the Poissonian noise. In all cases, the distribution is strongly biased 
towards the high-/? section, reaching a maximum aX (3 > 0.9. It is now clear that at the proto- 
halo sites the tidal shear and inertia momentum tensors are strongly correlated with each 
other, regardless of the smoothing scale , which confirms quantitatively the previous works 
dLee fc Penlboool : IPorciani et al1l2002bh . Note also that for the case of R^ = 0.5 h'^Mpc a 
small number of proto-halos exhibit exceptionally low values of (3 < 0.5, while for the cases 
of Rs = 1, 2, 5 h~^Mpc all proto-halos have f3 > 0.5. 

This result also implies that a halo of mass M tends to be made of those dark matter 
particles that accreted into the initial sites along the principal axes of the local tidal field 
smoothed on the same mass scale M, which is consistent with the Zel'dovich approximation 
dZerdovichlflQTol ). 



3.1. Dependence on Protohalo's Sphericity and Linear Density 

We investigate how (3 varies with the shapes of the proto- halos. Using the t hree eigen- 
values of I, we measure the sphericity S of each proto-halo as (IHahn et al.l 120071 1 



S 



1/2 



(3) 



where gi and represent the largest and the smallest eigenvalue of I, respectively. Binning 
the range of S, we measure the mean of /5 averaged over a given differential bin, [S, S + dS]. 
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Figure [2] shows the scatter-plot of (3 vs. S. The thick sohd hne corresponds to the mean value 
{(3) as a function of S. The errors represent one standard deviation in the measurement of {(3) 
calculated as ^J [n — 1) where n represents the number of the proto-halos belonging 

to a given differential bin [S*, S + dS] . It can be noted that (3 increases almost monotonically 
with 5*. In other words, the less spherical a proto-halo is, the weaker is the correlation 
between the tidal shear and the inertia momentum tensors. For those proto-halos with low 
sphericity {S < 0.4), the mean value, (/?), drops below 0.7. 

We also investigate how (3 changes with the local density field at the proto-halo sites. 
Let 5i denote the initial density contrast measured at the proto-halo's center of mass. Since 
the three simulations used here started at different initial redshifts Zi (see §2), we use the 
linearly extrapolated linear density 5 to 2; = instead of 5i itself, which is calculated as 
5l = [D{z) / D{{)y\5i where D{z) is the linear growth factor. We bin the range of 6 and 
calculate the mean of (3 averaged over a given differential bin, [5, 5 + d5\. Figure [3] shows the 
scatter-plots of (3 vs. 5. The thick solid line corresponds to the mean value (/?) as a function 
of 5. Note first that for the case of the large filtering radius Rs = 5h~^Mpc all values of S lie 
in quite a narrow raiige converging to the cri t ical density value 6p,r ~ 2.5 for the ellipsoidal 



collapse (jSheth et al.l 1200 ll : iDesjacquesI l2008l : iRobertson et all 120091 ). In contrast, for the 
cases of smaller filtering radii Rg = 0.5, 1, 2 h~^Mpc, the values of S spread over quite a 
wide range from —1 to 10. The value of f3 tends to be highest when 6 becomes close to 6ec- 
As 6 deviates from Sec, the value of (3 decreases both in the high-(5 and in the low-5 sections. 

Since S and 6 are not mutually independent but S tends decrease in the low-5 region, it 
may be interesting to see how (3 vary as 5* and 6 change simultaneously. Figure H] shows the 
contour-plots of P in the S-6 plane. This plot clearly demonstrates that those proto-halos 
with high S and 6 ~ 6ec tend to have a high value of /3. 



3.2. Dependence on Filtering Scale 

We now explore the dependence of /3 on the filtering radius. Let us consider two different 
scales Rth and Rg. We first calculate the inertia tensors of those halos whose masses lie in 
a fixed range of (0.9Mth, l.lMth) where Mth denotes the top-hat masses enclosed by the 
filtering scale -Rth- Then, we calculate the mean correlations of these inertia tensors with 
the tidal tensors smoothed on the different scale of Rs- For each given Rth, we consider four 
different scales Rg for the smoothing of the tidal tensors. 

Figure [5] shows as a function of the filtering scale Rg. Each panel plots the mean 
values of P calculated using the inertia tensors of the proto-halos with mass Mth enclosed by 
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a fixed radius Rth and the tidal tensors smoothed on four different scales Rs- It can be seen 
that increases monotonically when Rs < 2i?th and decreases sharply when Rg > 2Rt^. 
The value of (3 reaches a maximum when i?th < -Rs < 2i?th- In terms of mass, it can be said 
that f3 becomes maximal at Mth < Mg < SMth- In other words, the inertia tensors of the 
proto-halos with masses Mth are strongly correlated with the local tidal fields smoothed on 
the mass scales that have the same order of magnitude as Mth- If the tidal fields are smoothed 
on the scales order of magnitude smaller or larger than Mth, the correlations between the 
two tensors decrease. Thus, the correlation strength between the two tensors depends on 
the scale on which the tidal field is smoothed. 



IMPLICATION ON THE LOW SURFACE BRIGHTNESS GALAXIES 



In the light of our results, it is time to recall the linear tidal torque theory according to 
which the magnitude of the specific angular momentum (angular momentum per unit mass) 
of a proto -halo increases as the correlation of its inertia tensor with the local t i dal te nsor 
increases JPoroshkevichl ll97ol : Iwhitel llQsi ICatelan fc TheunsI Il99d : Ihee fc PenI boooh . In 
the context of the linear tidal torque theory, those proto-ha lo sites which have lo w values 
of (3 are likely to acquire higher specific angular momentum (jPorciani et al.ll2002al ). To test 
this core prediction of the linear tidal torque theory, we explore the dependence of /3 on the 
specific angular momentum measured in Eulerian space. For each halo, we first measure the 
angular momentum vector J in the Eulerian space as 



E 



(4) 



where and represent the position and the velocity of the a-th particle in the halo's 
center of mass frame, respectively. The specific angular momentum vector j of each halo is 
then calculated as j = J/M. Binning the range of j, we calculate the mean of jS averaged 
over a given differential bin, + dj]. Figure [6] shows the scatter-plot of (3 vs. j. The 
thick solid line corresponds to the mean value {I3){j) and the errors represent one standard 
deviation in the measurement of {l3){j). As can be seen, there is a strong trend that (3 
increases with j. In other words, those proto-halos whose inertia tensors are less correlated 
with the local tidal tensors are likely to acquire higher specific angular momentum vectors 
in the subsequent evolution. 

A crucial implication of our results is that the proto-halos with the lowest 5* will thus 
acquire the highest specific angular momentum since the proto-halos with the lowest S 
are found to have the weakest correlations between the inertia and tidal shear tensors (see 
Figure [2]). It is very interesting to recall that the high-specific angular momentum is the 
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characteristic property of the low-surface brightness galaxies (LSBGs) flBoissier et al.l 12003 



Monnier Ragaigne et al.ll2003l . and references therein). Our result leads to a speculation that 
the LSBGs might originate from those peculiar proto-halos with very low values of S whose 
inertia tensors are weakly correlated with the local tidal tensors. 

It is, however, worth mentioning here that the low values of 5* do not necessarily corre- 
spond to elongated particle distribution of proto-halo regions. A proto-halo region can have 
low values of 5* if it consists of two disconnected patches with center of mass located in the 
middle. To examine whether or not the proto-halos with low values of S have connected 
particle distribution, we inspect the particle distribution of those proto-halos with S below 
the tenth percentile in the mass range 2 < M / [10^^ Mq] < 4. We consider only the 
low-mass halos for this inspection since they exhibit the lowest values of S. It is found that 
among the inspected proto-halos, approximately two third of them have connected particle 
distribution while the other one third consist of more than two patches. 

Figure [7] shows an example of a connected proto-halo, plotting its particle distribution 
in the z-y , x-y and z-x plane in the top- left, top-right and bottom-left panel, respectively. 
The sphericity of this example is as low as S* = 0.24. As it can be seen, the Lagrangian region 
is connected, having indeed quite an elongated shape. Figure |8] shows an example of the 
disconnected proto-halo regions. Its sphericity is found to be S* = 0.25. Strictly speaking, 
the first order linear tidal torque theory is not valid to explain the generation of the angular 
momentum of such a disconnected proto-halo as shown in Figure [HI In this case it might 
be the gravitational merging of the disconnected patches in the subsequent evolution rather 
than the misalignments between the inertia and tidal tensors in Lagrangian space th at would 



contribute to the built-up of the higher angular momentum (jVitvitska et. al.ll2002l ). 



Whether or not a proto-halo region is connected, however, the particles of a proto-halo 
region with lower value of /3 (and thus lower values of S) will end up in a final halo with 
higher specific angular momentum, as revealed in Figure [61 Therefore, it is still possible 
to postulate that the peculiar proto-halos with low values of /5 correspond to the present 
LSBGs, no matter what caused the built-up of high angular momentum of the LSBGs. 



DISCUSSION AND CONCLUSION 



In the classical Zel'dovich model (IZerdovichlll970l ). the inertia tensors of bound objects 
are perfectly correlated with the local tidal tensors in Lagrangian space. In practice, the 



two tensors are found to be indeed strongly but not perfectly correlated (ILee fc Peril [2000 



Porciani et al.l l2002bl ). Here, we have determined quantitatively how the strengths of the 
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correlations between the two tensors depend on proto-halo's shape and mass, local density 
and filtering scale. Since the proto-halos form through the tidal flows of dark matter particles 
along the principal axes of the local tidal fields from the surrounding matter distribution, 
it is in fact natural to expect strong correlations between the two tensors. Deviations from 
the perfect correlations of the two tensors imply the existence of higher order perturbations 
from the simple tidal flows of CDM particles. 

We have also found that for the peculiar proto-halos with low-S" the correlations between 
the two tensors tend to be weak (i.e., having low value of P). Since those proto-halos which 
exhibit lower value of /? end up in halos with higher specific angular momentum, it is intrigu- 
ing to speculate a hypothesis that those peculiar proto-halos with lowest values of S would 
develop into the low surface brightness galaxies (LSBGs) at present epoch. Since the LSBGs 
are believed to be dark-matter dominated, t heir density profiles a re often directly compared 



with that of the dark halos (i.e., NFW profile iNavarro et al.lll996l ) and the shallow inner-core 



slope of the observed density prof iles of LSGB s has been used as a counter-evidence f or the 



cold dark matter paradigm (e.g., lMoorelll994l ). For instance, iKuzio de Naray et al.l (120081 ) 



have recently studied the rotation curves of 17 LSBGs obtained from the high-resolution op- 
tical velocity fields from DensePak spectroscopic observations and shown that the observed 
LSBGs can be matched with the NFW halos only if the LSBGs have 20 kms~^ non-circular 
motions. If the LSBGs originate from the peculiar proto-halo sites with exceptionally low 
value of (3 as in our hypothesis, then their proto-halo sites may have had very low sphericity 
as our results imply. Those peculiar proto-halo sites which have extremely low sphericity at 
the initial stages might as well develop non-circular motions. 

It will be interesting to study numerically the density profiles of those proto-halos with 
exceptionally low values of (3 and compare them with the standard NFW ones. It will be also 
interesting to compare the number density of those proto-halos with low values of f3 with 
that of the observed LSBGs as a function of mass. We plan to work on these two projects 
and hope to report the results elsewhere in the future. 
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Fig. 1. — Probability density distribution of (3 of the proto-galactic sites with Poissonian 
errors when the initial density field is smoothed on the scale of 0.5, 1, 2 and 5/i~^Mpc 
(top-left, top-right, bottom-left, and bottom-right, respectively). 
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Fig. 2. — Scatter plots of (3 vs. the sphericity S of the proto-halo sites when the initial 
density field is smoothed on the scale of 0.5, 1, 2 and 5 h~^Mpc (top-left, top-right, bottom- 
left, and bottom-right, respectively). In each panel, the solid line represents the mean values 
and the errorbars indicate the standard deviation in the measurement of (jS). 
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Fig. 3. — Scatter plots of jS vs. the linearly extrapolated density 6 of the proto-halo sites 
when the initial density field is smoothed on the scale of 0.5, 1, 2 and 5 /i~^Mpc (top- 
left, top-right, bottom-left, and bottom-right, respectively). In each panel, the solid line 
represents the mean values and the errors correspond to the one standard deviation in 
the measurement of 
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Fig. 4. — Contour plots of (3 in the S-5 when the initial density field is smoothed on the scale 
of 0.5, 1, 2 and 5/;,~^Mpc (top-left, top-right, bottom-left, and bottom-right, respectively). 
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Fig. 5. — Mean values as a function of the filtering radius Rg for four different values 
of the proto-halo masses. The four masses are enclosed by the four top-hat radii i?th as 
Mth = (47r/3)p-Rj^jj where i?th = 0.5, 1, 2 and 5 /i~^Mpc (top-left, top-right, bottom-left, and 
bottom-right, respectively). 
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Fig. 6. — Scatter plots of /5 vs. the specific angular momentum j of the proto-halo sites 
when the initial density field is smoothed on the scale of 0.5, 1, 2 and 5h~^Mpc (top- 
left, top-right, bottom-left, and bottom-right, respectively). In each panel, the solid line 
represents the mean values and the errors correspond to the one standard deviation 

in the measurement of {(3){j). 
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Fig. 7. — Particle distribution of a proto-halo with S = 0.24 in the x-y , y-z and x-z plane 
in the top- left, top-right and bottom-left panel, respectively. 
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Fig. 8. — Same as Figure [7] but for the case of a disconnected proto-halo with 5" = 0.25. 
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Table 1. Filtering radius (Rs), the top-hat mass (M^) enclosed by Rg, the number of halos 

(Nh) with mass in range of [0.9Ms, l.lMs] 



Rs 




Nh 








0.5 


3.6 X 10^° 


1049 


1 


2.9 X 10" 


1472 


2 


2.3 X 10^2 


1835 


5 


3.6 X 10^3 


145 



